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Abstract 
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(N 

HH Cosniogenic neutrinos originate from photo-hadronic interactions of cosmic ray protons with the cosmic 

ri microwave background (CMB). The neutrino production rate can be constrained through the accompanying 

Q^ electrons, positrons and 7-rays that quickly cascade on the CMB and intergalactic magnetic fields to lower 

JL energies and generate a 7-ray background in the GeV-TeV region. Bethe-Heitler pair production by protons 

\^ also contributes to the cascade and can tighten the neutrino constraints in models where extragalactic cosmic 

j/j rays begin to dominate over the galactic component at a relatively low "crossover" energy. We investigate 

Ci this issue in the light of the recent Fcrmi-LAT measurements of the diffuse extragalactic 7-ray background 
and illustrate by a fit to the HiRes spectrum how the prediction of the cosmogenic neutrino flux in all-proton 

^SJ models varies with the crossover energy. The neutrino flux is required to be smaller when the 7-ray bound is 

^ applied, nevertheless such models are still consistent with HiRes and Fcrmi-LAT if one properly takes into 

^—^ account the energy uncertainty of cosmic ray measurements. The presently allowed flux is within reach of 

S^ the IceCube neutrino telescope and other dedicated radio experiments. 
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1. Introduction 

Soon after the discovery of the cosmic microwave background (CMB) [T], it was reahzed that interactions 
of extragalactic ultrahigh energy (UHE) cosmic rays (CRs) on the rehc photons would suppress the cosmic 
ray flux at energies > 5 x 10^° GeV, the so-called "GZK cutoff" [2 [3]. It was pointed out subsequently [3] 
that the GZK interaction also generates a "cosmogenic flux" of neutrinos, through the decay of secondary 
charged pions. Forty years later, the predicted suppression of the UHE CR flux was indeed observed by the 
HiRes [S] and Auger [|6[ experiments. However, the cosmogenic flux of neutrinos has yet to be detected. 

The GZK reaction chain generating cosmogenic neutrinos is well known [7]. The intermediate state of 
the reaction P7cmb — >■ mr'^ /pn^ is dominated by the A"*" resonance, because the neutron decay length is 
smaller than the nucleon mean free path on the CMB. Resonant pj interactions produce twice as many 
neutral pions as charged pions. Direct pion production via virtual meson exchange contributes only about 
20% to the total cross-section, but is almost exclusively into 7r+. Hence, p'y interactions produce roughly 
equal number of tt"*" and 7r°. Gamma-rays, produced via 7r° decay, subsequently cascade electromagnetically 
on intergalactic radiation fields through e~^e~ pair production followed by inverse Compton scattering. The 
net result is a pile up of 7-rays at GeV-TeV energies, just below the threshold for further pair production 
on the diffuse optical background. Meanwhile each 7r+ decays to 3 neutrinos and a positron; the e+ readily 
loses its energy through inverse Compton scattering on the diffuse radio background or through synchrotron 
radiation in intergalactic magnetic fields. The neutrinos carry away about 3/4 of the tt"*" energy, therefore 
the energy in cosmogenic neutrinos is about 3/4 of that produced in 7-rays. 

The normalization of the neutrino flux depends critically on the cosmological evolution of the CR sources 
and on their proton injection spectra [8j. It also depends on the assumed spatial distribution of sources; for 
example, local sources in the Virgo cluster [5], would dominate the high energy tail of the proton spectrum. 
Another source of uncertainty is the energy at which there is a transition from Galactic to extragalactic 
CRs as inferred from a change in the spectral slope. The "ankle" at ~ 3 x 10^ GeV seems to be a natural 
candidate for this transition ITUIIIIIIII], but a lower energy crossover at the "second knee" at ^ 5 x 10^ GeV 
has also been advocated [13l [14] . A fourth source of uncertainty is the chemical composition of the parent 
CRs - if these are heavy nuclei rather than protons, then the neutrino flux is reduced ^5j. 

The most up-to-date calculation |16j of the cosmogenic neutrino flux combines a double-fit analysis of 
the energy 117] and elongation rate |18| measurements to constrain the spectrum and chemical composition 
of UHE CRs at their sources. Injection models with a wide range of chemical compositions are found to 
be consistent with observations. The data is consistent with a proton-dominated spectrum with a small 
admixture of heavy nuclei, in which case the cosmogenic neutrino flux is rather similar to the all-proton 
model. In this case, kilometer-scale neutrino telescopes are expected to observe of 0(1) cosmogenic neutrino 
event per year. In contrast, an intermediate to heavy nuclear composition beyond the ankle, as indicated by 
the elongation rate vs. energy [TBI, can lead to a considerable suppression (up to two orders of magnitude) 
of the cosmogenic neutrino flux in comparison to the all-proton case. However, little is known about the 
chemical composition from just below to beyond the GZK cutoff, where the most significant contribution to 
cosmogenic neutrinos form UHE CR protons is expected. It is interesting to note that uncertainties in the 
extrapolation of the proton-air interaction - cross-section, elasticity and multiplicity of secondaries - from 
accelerator measurements to the high energies characteristic for air showers are large enough to undermine 
any definite conclusion on the chemical composition (THl HO] . 

In this work we study the constraint set by the diffuse 7-ray background on all-proton models of extra- 
galactic CRs. We parametrize our ignorance of the crossover energy - which marks the transition between 
the galactic and extragalactic components - as a variable low energy cutoff in the proton injection rate. By 
fitting only to CR data above the crossover energy, taken to be between 10^^'^ eV and 10^^ eV, we determine 
the statistically preferred values of the spectral index 7 and cosmic source density evolution index n hy a 
goodness-of-fit (GOF) test of the HiRes data, taking into account the energy resolution of about 25%. For 
each model we check that the total energy density of the EM cascade is below a critical value inferred from 
the recent measurement of the extragalactic 7-ray background by the Fermi-LAT Collaboration [5T] . We 



find that the allowed range of the cosmogenic neutrino flux increases with the crossover energy and can be 
up to an order of magnitude larger than the values presented in a recent study [2 2) . 

This paper is organized as follows. We begin in §[2] with a discussion of the extragalactic proton fluxe 
and the corresponding energy density of cascade 7-rays. We present our statistical method in §[3] and discuss 
our results in § |4] 

2. Extra-galactic Proton Fluxes and Diffuse Gamma Background 

For a spatially homogeneous distribution of cosmic sources, emitting UHE particles of type «, the co- 
moving number density Yi is governed by a set of (Boltzmann) continuity equations of the form: 

Y, = dE{HEYi) + dEihYi) -r,Y, + J2j d^j 7r<^j + ^^ : (1) 

i 

together with the Friedman-Lemaitre equations describing the cosmic expansion rate H{z) as a function of 
the redshift zjj The first and second terms on the r.h.s. describe, respectively, redshift and other continuous 
energy losses (CEL) with rate b = dE/dt. The third and fourth terms describe more general interactions 
involving particle losses (i — >■ anything) with interaction rate F^, and particle generation of the form j ^ i 
with differential interaction rate jij. The last term on the r.h.s., Ct, corresponds to the emission rate per 
CO- moving volume of CRs i. We refer to Ref. [21] for explicit definitions of the coefficients in Eq. (IT]). 

Extragalactic protons lose their energy via Bethe-Heitler (BH) pair production and photo-hadronic in- 
teractions on cosmic radiation backgrounds, notably the CMB. Bethe-Heitler pair production, p + 7bgr — >■ 
p+e++e^, can be treated as a continuous energy loss due to its low inelasticity [25]. This process dominates 
the evolution of the spectra at energies between 10^ GeV and a few times 10^° GeV. At higher energies, 
resonant photo-hadronic interactions with CMB photons lead to a sharp suppression of the spectrum [U [3] . 
The produced charged and neutral pions release electrons, positrons, neutrinos and photons through their 
decay. We calculate the spectra of hadrons and neutrinos using the Monte Carlo package SOPHIA [^ . 

It is possible to approximate the energy loss in the hadronic cascade due to photo-pion production as a 
CEL with 

^(z,i?) = b{z,E) ~ ETp{z,E) - JdE'E'jppiz,E,E'). (2) 

Diffractive pj processes at high energies with large final state multiplicities of neutrons and protons ulti- 
mately invalidate the CEL approximation. However, the relative error below 10^^ GeV is less than 15% so 
we will use this approximation for a detailed numerical scan in the model space of proton spectra. 

Electromagnetic (EM) interactions of photons and leptons with the extragalactic radiation backgrounds 
and magnetic field can happen on time-scales much shorter than their production rates. The relevant 
processes with background photons contributing to the differential interaction rates 7ee, J-ye and j^^y are 
inverse Compton scattering (ICS), e^ + 7bgr — > e* + 7, pair production (PP), 7 -I- 7bgr — >■ e+ + e^, double 
pair production (DPP) 7 -I- 7bgr — ?> e+ + e^ + e+ + e^, and triple pair production (TPP), e* + 7bgr — > 
e^ -f e+ + e~ [25l [27] . High energy electrons and positrons can also lose energy via synchrotron radiation on 
the intergalactic magnetic fieldrlthe strength of which is limited to be below ~ 10~^G [30] and suggested to 
be of O(10^^^)G by simulations of large-scale structure formation [32j; details on the calculation are given 
in [Appendix A] 

^This is given by H^{z) = Hq [f2ni(l + z)^ + S^a], normalised to its value today of Hq ~ 70 kms"^ Mpc~^, in the usual 
"concordance model" dominated by a cosmological constant with Q\ ^ 0.7 and a (cold) matter component, Qm ~ 0.3 I23| . 
The time-dependence of the redshift can be expressed via dz = — dt (1 -|- z)H. 

^ Since we consider a relatively strong intergalactic magnetic field, we can neglect TPP by electrons in the following I28| . 
Also below lO^'^ GeV we can safely neglect DPP of photons in the calculation |29| . 



Synchrotron radiation in strong magnetic fields can also by-pass the EM cascade and transfer energy to 
sub-GeV photons that are unconstrained by the Fermi-LAT spectrum [31]. In the case of a strong 10~^G 
field this can be relevant for electrons around 10^ GeV, where synchrotron loss starts to dominate over 
ICS loss in the CMB and where the corresponding synchrotron spectrum still peaks below 100 MeV. In 
our calculation we adopt a moderate value of 10~^^G following [35]. However, we have checked that a 
significantly larger field strength of lO^^G has little effect on the 7-ray flux in the Fermi-LAT energy range 



relevant for our discussion (see Fig. C.8 in the Appendix). 

In this paper the emission rate of CR protons per co-moving volume (GeV~^ cm~'^ s~^) is assumed, as 
per usual practice, to follow a power-law: 

{/_(£'/£'i„in) E < -Bmin , 
1 E^i^ <E< ^^ax , (3) 

/+(-^/-^max) £'max<-E. 

We will consider spectral indices 7 in the range 2-^3. The functions f±{x) = x^^ exp(l — x^"^) in Eq. (|3| 
smoothly turn off the contribution below E'min and above -Emax- We set E'max = 10^^ eV in the following 
and vary -Emin in the range 10^^'^ -;- 10^^ eV, corresponding to a galactic-extragalactic crossover between the 
"second knee" and the "ankle" in the CR spectrum. 

The cosmic evolution of the spectral emission rate per comoving volume is parameterized as: 

Cpiz,E)=n{z)Cp{0,E). (4) 

For simplicity, we use the standard approximation 

n{z) = {l + zrQ{z„,,^-z), (5) 

with 2;max — 2. Note that the dilution of the source density due to the Hubble expansion is taken care of 
since C is the comoving density, i.e. for no evolution we would simply have Ji ^ I. We consider cosmic 
evolution of UHE CR sources with n in the range 2 -;- 6. 

As mentioned above, EM interactions of photons and leptons with the extra-galactic background light 
and magnetic field can happen on time-scales much shorter than their production rates. It is convenient to 
account for these contributions during the proton propagation as fast developing electro-magnetic cascades 
at a fixed redshift. We will use the efficient method of "matrix doubling" [33] for the calculation of the 
cascades. Since the cascade 7-ray flux is mainly in the GeV-TeV region and has an almost universal shape 
here, it is numerically much more efficient to calculate the total energy density Wcas injected into the cascade 
and compare this value to the limit imposed by Fermi-LAT. The total energy density (eV cm~^) of EM 
radiation from proton propagation in the past is given as 

C^eas = J dEEncas(0, E) ^ J dt J dE ^J^^^ ""P^^' ^^ ' ^^^ 

where n{z, E) is the physical energy density at redshift z, defined via n{z, E) = {1 + z)^Y{z, E). We discuss 
the derivation of this equation in Appendix B The continuous energy loss of protons into the cascade. 



denoted by 6cas7 is in the form of electron, positron and 7-ray production in BH (^bh) and photo-pion (foir) 
interactions 13 

In the following we derive the BH and photo-pion contribution to ojcas separately. For the photo-pion 
contribution we estimate 

b^{z, E)^ j dE'E' [7pe- (z, E, E') + 7p,+ (z, E, E') + 7p^(z, E, E')] (7) 



^Note the difference between our Eq. ([ml and Eq. (10) in Ref. 122' where the approximation dEb{z, E) ~ b{z, E)/ E is used 
(c/. Fig. 5 in Ref. [24]) and an adiabatic scaling with redshift, b{z, £) ~ (1 -(- 2)^6(0, E(l + 2)), is assumed. 
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Figure 1: Left Panel: Goodness of fit test of the HiRes data [5]. We show the 68% (pink), 95% (blue) and 99% (magenta) 
confidence levels of the injection index 7 and the cosmic evolution index n. The black lines indicate the allowed regions before 
the cascade (ojcas) bound is imposed. Right Panel: The corresponding energy density in the EM cascade. 

For the energy loss via BH pair production we use the expression given in Ref. [53]. Note, that since the 
photo-pion contribution in the cascade is dominated at the GZK cutoff its contribution should increase 
relative to BH pair production with increasing crossover energy and, hence, also the associated neutrino 
fluxes after normalization to 7-ray and CR data. 



3. Goodness of Fit Test 



In this section we present the results of a GOF test of the compatibility of a given model, characterized by 
the injection spectral index 7, cosmic evolution index n, and crossover energy -Emin, with the CR experimental 
data - in particular HiRes I and H - imposing also consistency with the Fermi-LAT measurements of the 
diffuse extra-galactic 7-ray background. 

Given the acceptance A^ (in units of area per unit time per unit solid angle) of the experiment for the 
energy bin i centered at Ei and with bin width A^, and the energy scale uncertainty of the experiment, aE^ 
the number of expected events in the bin is given by 



E,{l+5)+Ai/2 



m{n,^,M,S)=A, 



JknniE)dE, 



(8) 



Ei{l+S)-Ai/2 



where J^ „ ^{E) = np{0, E)-^ is the proton flux arriving at the detector corresponding to a proton source 
luminosity as in Eq.Jsl), with the cosmic evolution of the source density given by Eqs. ^ and (Is]). The 
parameter 5 in Eq. (|8| above is a fractional energy-scale shift that reflects the energy-scale uncertainty of 
the experiment, and A/" is the normalization of the proton source luminosity. 

The probability distribution of events in the z-th bin is of the Poisson form with mean Ni. Corre- 
spondingly the r-dimensional (r being the number of bins of the experiment with Ei > -Emin) probability 
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Figure 2: The allowed proton flux (at the 99% confidence level) for increasing crossover energy -Eniin- Each fit of the proton 
spectrum is marginalized with respect to the experimental energy uncertainty and we show the shifted predictions in comparison 
to the HiRos central values [5]. For comparison we also show the Auger data [6l ll7| which has not been included in the fit. 



distribution for a set of non-negative integer numbers k = {fci, ...k,.}, Pt{n,"f,Af,S), is just the product of 
the individual Poisson distributions. 

According to this r-dimensional probabiHty distribution, the experimental result iV™P = {A''™'', ..., N^^^} 
has a probability -Pj^cxp (n, "f,J^, S) and correspondingly the experimental probability after marginalizing over 
the energy scale uncertainty and normalization is: 



^oxp (?^, 7) = Max5,AA^«oxp (n, 7, A^, (5) 



(9) 



where the maximization is made within some prior for 6 and Af. For the energy shift S we have used two 
forms for the prior, either a top hat spanning the energy-scale uncertainty of the experiment, aE^ , or a 
gaussian prior of width a^s ■ 

For J\f wc impose the prior arising from requiring consistency with the Fermi-LAT measurements |21) 
of the diffuse extra-galactic 7-ray background. In order to do so we obtain the total energy density of EM 
radiation from the proton propagation using Eq. ([6]) and we require following Ref. |22| : 

WcasCA/", n, 7) < 5.8 X 10"^ eV/cm^ . (10) 

The marginalization in Eq. ([9| also determines Abest a-nd (5bcst for the model, which are the values of the 
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Figure 3: Systematic effect of the experimental energy resolution on the fitted spectral index 7 and cosmological evolution 
parameter n. For illustration we show the dependence of the 95% C.L. bound for a crossover energy of 10^* eV. The blue 
contour corresponds to the region shown in Fig. IT] assuming an uncorrelated energy shift of 25% in both data sets (HiRes I 
and II) \5\ , for a flat prior ( "top-hat" distribution)"^ The red dashed curve assumes correlated errors of the energy resolution in 
both data sets. The black dotted curve shows the result for uncorrelated errors with a Gaussian prior, and the dashed-dotted 
line shows uncorrelated errors with a flat prior, but with a lower uncertainty of 15%. 

energy shift and normalization that yield the best description of the experimental CR data, subject to the 
constraint imposed by the Fermi-LAT measurement. 

Altogether the model is compatible with the experimental results at given goodness of the fit (GOF) if 
^F^(n,7,A/'best,'^bcst)e [Pj:(n,7,A/'best,'5bcst) - ^oxp(?^,7)] < GOF (11) 



Technically, this is computed by generating a large number A'lop of replica experiments according to the prob- 
ability distribution Pt{n, 7, A/bosti <5bcst) and imposing the fraction F of those which satisfy Pt{f^i It^i ^bcst) > 
-Pexp(",7) tobeF< GOF. 

With this method we determine the value of (n, 7) parameters that are compatible with the HiRes I and 
HiRes H experiments [5]. We plot in the left panel of Fig. IT] the regions with GOF 64%, 95% and 99% for 
four values of the minimum (i.e. crossover) energy. In the right panel we show the corresponding ranges 
of Wcas,best for thc modcls as a function of the cosmic evolution index n. In order to display explicitly the 
impact of the constraint from the Fermi-LAT measurements of the diffuse extra-galactic 7-ray background 
(10 1, we show the corresponding GOF regions without imposing that constraint. In Table. Il| we list the 
parameters corresponding to the best-fit models and to the models with minimal and maximal contributions 
to Wtt and Wcas ~ i-^tt + i^BH at the 99% C.L., together with the corresponding energy shifts which give best 
fits to the HiRes I and Hires II data. We also show the parameters for the models with maximum Wtt and 
Wcas without imposition of the Fermi-LAT constraint. 

As an illustration of thc agreement with the CR data we show in Fig. [T] the range of proton fluxes 
corresponding to models with GOF 99% or better for increasing crossover energies i^min- As discussed above 
each fit of the proton spectra is marginalized with respect to the experimental energy scale uncertainty and 
we show thc shifted predictions with (Jbost hi comparison to the HiRes data at central value. We also show 
in the figure the results from Auger [SI HZ], though these have not been included in the analysis (see below). 
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Table 1: Cosmic ray source parameters which best fit the HiRes data [5], along with those which yield minimal and maximal 
contributions to uj-,^ (i.e. neutrino fluxes) and oJcas = tJvr -I- tJBH (i-O- 7-ray fluxes), all at the 99% C.L. 



These results are obtained assuming an energy scale uncertainty ue, = 5% with a "top-hat" prior for 
the corresponding energy shifts which are taken to be uncorrelated for HiRes I and HiRcs II. In Fig. [3] we 
explore the dependence of the results on these assumptions by using a different form for the prior, assuming 
the energy shifts to be correlated between the two experiments, or reducing the uncertainty to ue^ — 15%. 
As seen in the figure, the main effect is associated with the reduction of the energy scale uncertainty which, 
as expected, results in a worsening of the GOF for models with larger n. This is directly related to the 



normalization constraint from Eq. (10 1. If one naively ignores the energy scale uncertainty, the constraint 
in Eq. (10) rules out models with n ^ 3 (the precise value depending on the assumed -Emin)- However, once 



the energy scale uncertainty is included, the constraint of Eq. ( 10 ) plays a weaker role on the determination 
of the GOF of the models. It does however imply a maximum value of A/bcst which, as we will see, impacts 
the corresponding ranges of neutrino fluxes. 

The corresponding range of 7-ray and cosmogenic neutrino fluxes (summed over flavour) is shown in 
Fig. [4] for models with minimal and maximal energy density at the 99% C.L. As expected, the maximum 
7-ray fluxes are consistent with the Fermi-LAT data within the errors. For illustration, we also show as a 
dotted line the "naive" 7-ray limit E^Jcas ^ ccj™f'^/47rlog(TeV/GeV), corresponding to a 7-ray flux in the 
GeV-TeV range which saturates the energy density (10). 



We have not included in the analysis the results from the Auger Collaboration [SI [T7] , which are shown 
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Figure 4: Comparison of proton, neutrino and 7-ray fluxes for different crossover energies. We show the best-fit values (solid 
lines) as well as neutrino and 7-ray fluxes within the 99% C.L. with minimal and maximal energy density (dashed lines). The 
values of the corresponding model parameters can be found in Table. IT] The dotted line labeled "maximal cascade" indicates 
the approximat e lim it E^Jczs < ca;™|''/47rlog(TeV/GeV), corresponding to a 7-ray flux in the GeV-TeV range saturating the 
energy density l |10| . The 7-ray fluxes are marginally consistent at the 99% C.L. with the highest energy measurements by 
Fermi-LAT. The contribution around 100 GeV is somewhat uncertain due to uncertainties in the cosmic infrared background. 

in Fig J2] for illustration only (hence our results are directly comparable to those in Ref.[22|). As described 
in Refs. [HIIIZIj besides the energy scale uncertainty there is also an (energy-dependent) energy resolution 
uncertainty which implies that bin-to-bin migrations influence the reconstruction of the flux and spectral 
shape. Since the form of the corresponding error matrix is not public, this data E] cannot be analysed 
outside the Auger Collaboration. 



4. Discussion 



The cosmogenic neutrino fluxes that we have shown in Fig. [4] are compared to present upper limits on 
the diffuse neutrino flux in Fig. p] As before, the solid green line shows the neutrino flux (summed over 
flavours) corresponding to the best fit of the proton spectra and the dashed green line indicate the range of 
neutrino fluxes within the 99% C.L. For all crossover energies considered, the range of models at the 99% 
C.L. is consistent with existing neutrino limits. For illustration, the thin dotted line shows the larger range 
of neutrino fluxes at the 99% C.L. corresponding to a fit without the Fermi LAT constraint (c/. the black 
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Figure 5: The predicted best fit (solid) and 99% C.L. range of cosinogenic neutrino fluxes with (dashed) and without (dotted) 
the Fermi-LAT constraint. The values of the corresponding model parameters can be found in Table. IT] For comparison 
we show upper limits on the total diffuse neutrino flux from AMANDA 1351 I36| . Auger |34l . Lake Baikal 1371 . HiRes |38| 
(minimum of i/^ and u^ channel), RICE |39| and ANITA I40| . The black solid line shows the Sc sensitivity of IceCube after 
just 1 year of observation |42l . The cutoff at 10* GeV is artificial so we also show an extrapolation to higher energies as 
a black dashed line following Ref. l4in . AH limits are obtained assuming an equal distribution between neutrino flavours: 
A^i/j : N^^ : N^^ ~ 1 : 1 : 1 (and scaled appropriately where necessary). Integrated limits assuming an E~^ spectrum are 
shown as solid lines and differential limits as dotted lines (both limits are shown for Auger and ANITA). 

contours in the left panel of Fig. II]). It is apparent that this indirect bound from GeV-TeV 7-rays does 
reduce the number of possible models significantly. 

At this point it is worth stressing that the Fermi-LAT spectrum used in this analysis is not the result of a 
direct observation but is derived by a foreground subtraction scheme. The extra-galactic 7-ray background 
inferred by EGRET [H] shows a significantly larger intensity and a harder spectral index. A possible 
source of the differences could be due to the different diffuse galactic emission (DOE) models used in the 
analysis. As pointed out in \2\ a re-analysis of the EGRET data with an updated DGE model [45] is 
comparable with the intensity observed with Fermi-LAT. It is beyond the scope of this paper to address 
these systematic uncertainties. The maximal effect of a larger 7-ray background intensity is indicated by the 
extended parameter regions shown in Fig. [T] which are derived without the Fermi-LAT constraint together 
with the corresponding range of neutrino fluxes in Fig. [5] 

The overall range of neutrino fluxes increases along with the crossover energy - not only in magnitude, 
which is expected already due to the reduced set of CR data used in the GOF test, but also to significantly 
larger neutrino fluxes. Also the cosmogenic neutrino flux of the best-fit models increases by over a factor 
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of two in the peak region (~ 10^ GeV). This confirms our earher suspicion that an increasing value of the 
crossover energy aUows a larger contribution of cosmogenic neutrinos relative to the 7-rays and hence larger 
neutrino fluxes. 

Figure p] also shows the estimated sensitivity of IceCube |45j (5cr) to neutrino fluxes in the 10^-10* GeV 
(solid) |42)~and the lO^-lO^*^ GeV (dotted) [JTj energy range after one year of observation. IceCube located 
at the South Pole is presently the largest neutrino telescope. On completion in early 2011 it will consist 
of a km'^-scale detector of transparent glacial ice, that is constantly monitored for Cerenkov light emission 
of secondary charged particles from high energy neutrino interactions. It is apparent from Fig. [5] that 
IceCube's sensitivity after one year is already sufficient to probe cosmogenic neutrino fiuxes from an all- 
proton spectrum of extra-galactic cosmic rays. If the crossover energy exceeds lO"'^^''^ eV, the best-fit model 
of the HiRes data is within reach of IceCube. 

In summary we find that while the expected range of cosmogenic neutrino fluxes in all-proton models 
is indeed reduced due to the constraint from the diffuse 7-ray background measurements, neutrino fluxes 
compatible at 99% C.L with the HiRes and Fermi-LAT results can be larger than those presented in Ref.p2] 
by up to factor of ^ 30 for the same values of E'max = 10^^ eV and Zmax — 2. In particular the allowed 
cosmogenic flux is still within reach of neutrino observatories like IceCube. Furthermore, our results are 
obtained with the simple parametrization of the source spectral emission rate in Eq. ([3]); larger neutrino 
fluxes might be allowed with a more general spectrum than a simple power-law. 

One can also turn this argument around and use observation or non-observation of cosmogenic neutrinos 
in the near future to provide additional constraints on the composition of cosmic rays [24" . We have assumed 
here an all-proton composition for extra-galactic cosmic rays. However, as mentioned already, the chemical 
composition of UHE CRs is rather uncertain and may well be dominated by heavy nuclei. In this case the 
limits on diffuse neutrino fluxes can still serve as a probe of the possible proton fraction in cosmic rays |24l 146) 
and the limits from diffuse 7-rays serve as an additional probe |47j . 
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Appendix A. Cascade Solution 

The relevant processes with background photons contributing to the differential interaction rates "fee-, 
7.ye and je-i are inverse Compton scattering (ICS), e* + 7bgr — >■ e* -I- 7, pair production (PP), 7 -|- 7bgr — >■ 
e+ + e^, double pair production (DPP) 7 -f- 7bgr -> e^" + e^ + e+ + e~, and triple pair production (TPP), 
e^ +7bgr — > 6=*= +e+ +e~ [551 [27] • The angular-averaged (differential) interaction rate, P,; (7^^) is defined as 

1 

T^{z,E,)^- d cos 9 de{l-P cos 0)n^{z,e)al°\ (A.l) 

-1 

j,,{z,E,,E,) = r.,{z,Ei) ^{E,,E,) , (A.2) 
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Figure A. 6: The energy spectrum of the CMB |23| and the CIB in the IR/optial .48. and radio .54, range at z = 0. The thin 
dashed Hne shows our extrapolation to UV energies. 

where n^(z, e) is the energy distribution of background photons at redshift z and dNij/dEj is the angular- 
averaged distribution of particles j after interaction of particle i. 

Besides the contribution of the CMB the shape of the cascade spectrum depends on the cosmic in- 
frared/optical background. We will use a recent estimate ^3S] (that we extrapolate slightly to UV energies 



as seen in Fig. A. 6 1 and assume a redshift dependence following the star formation rate as described in 
Ref. [23]. This is consistent with the constraints on the 7-ray opacity of the universe set by HESS [i^ . 
MAGIC [50] and Fermi-LAT [5l]. 

We have little direct knowledge of the cosmic radio background. An estimate made using the RAE 
satellite [S3] is often used to calculate the cascading of UHE photons [551. A theoretical estimate has 
been made |54| of the intensity down to kHz frequencies, based on the observed luminosity function and 
radio spectra of normal galaxies and radio galaxies although there are large uncertainties in the assumed 
evolution. The calculated values are about a factor of ~ 2 above the measurements and to ensure maximal 
energy transfer in the cascade we will adopt this estimate and assume the same redshift scaling as the cosmic 
infrared/optical background. The magnitude of the adopted radio background radio is not important for 
the shape of the GeV-TeV spectrum as can be seen from Fig. 6 of Ref. [55^ where even higher values are 
considered. We summarize the adopted cosmic radiation backgrounds in Fig. |A.6 



High energetic electrons and positrons may also lose energy via synchrotron radiation in the intergalactic 
magnetic field B with a random orientation sin 9 with respect to the velocity vector. We will assume in the 
following that the field strength B is of 0(10-^^)0 [32J, which leads to an efficient transfer of energy into 
the EM cascade. The synchrotron power spectrum (W eV~^) has the form 



_ , ^ v3a eB sin 9 , , , , , , 

V{E,,E^) = ^ -—F{E^/E,); F{t) = t / dz 1^5/3 W, 



(A.3) 



where we follow the notation of Ref. [56] with Ec = {3eB sin9 /2me){Ee/me)'^ . This can be treated as a 
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continuous energy loss of the electrons and positrons with a paranieteir] 

bsyn{Ee) = ^ / d COS fdE^ViE^E^) 



4a feBV f E,^ "^ 



9 V m, 



(A.4) 



We will assume in the calculation that the intergalactic magnetic field is primordial with a (fiux-conserving) 
redshift dependence B{z) = (1 + z)'^B{{)). Note, that the synchrotron energy loss has then a redshift 
dependence similar to BH pair production in the CMB, i.e. bsyn{z,E) = (1 + z)^6syn(0, (1 + z)E). It 
is also convenient to define jl^'^{Ee,Ej) = V{Ee,E^)/Ej, which has an analogous redshift dependence 
I.e. 7|y"(z, E,, E^) = (1 + zY^f^^{Q, (1 + z)E„ (1 + z)E^). 

The fast evolution of the cascade is governed by the set of differential equations, 

diY^{E) =-T^{E)Y-,{E)+ j AE' '^^^'^^\ ,{E')+ j dE'^,^{E',E)Y,{E') , (A.5) 

diY,{E) = -T,{E)Y,{E) + dE{b{E)Y,{E)) + j &E' Yl.,,{E' ,E)Y^{E') + ^,,{E' ,E)Y,{E')] , (A.6) 

which determines the evolution on s/iort time-scales AiFp^ <C 1 (the redshift z is kept /lied meanwhile). The 
initial condition Y^ / ^{E)\^^q is given by the sum of previously developed cascades and the newly generated 
contributions from proton interactions. 



The solution of Eqs. (A.5) and (A.6) for an infinitesimally small step Ai can be written for a discrete 



energy spectrum, Ni ~ AEiYi, as 



(UAh-S^f^-^^^^^^ re,(AiA 



r^e(Af) T,,{Ai)J^\N, 



N^ 



(i). 



(A.7) 



With the transition matrix T(Ai), defined by Eq. (A.7|, we can efficiently follow the development of the 
EM cascade over a distance At = 2^ At via matrix doubling [33 : 



r(2~A£) ~ [r(Ai)] 



N+l 



(A.; 



We will compare our calculation with results from other investigations in Appendix [Appendix C[ 



Appendix B. Energy Density of the Cascade 

We can express the system of partial integro-differential equations (fTl) as, 

Z, = ds {h{z,£)Z,{z, E)) - T,{z,£) Z, + {1 + z)Cf{z,£) 



(B.l) 



where we have defined £ — [1 + z)E, and Zi(z,E) = (1 + z)Yi{z,£), subject to the boundary condition 
Zj{Zma.x,E) = 0. The effective source term in Eq. ( |B.1[ ) is 

Cf ^Q + Y, f dE, 7,, (z, £j , E,) Zj , (B.2) 

The total energy of the cascade is given in Eq. ^. This can be obtained by integrating Eq. (B.l): 

JdEEZ,,,iz,E) =-JdEEd£[b,,,iz,£)Zp{z,E)]. (B.3) 



d 
dt 



*Note, the identity J dE [E dEibn^) + J dE''P{E' , -E)ne] = 0, implying overall energy conservation. 
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Figure B.7: Comparison of our calculations with the 7-ray spectra of Ref. |57| shown as black doted lines. We normalize the 
7-ray (green) and electron/positron (red) spectra to the proton spectra (black solid line). For comparison, the dashed lines 
show the spectra without the contribution of BH pairs and the dashed-dotted lines show the spectra without EM cascades. 



Integrating the r.h.s. by parts yields 



r.h.s. 



d^^F 



E- 



l + z 



bci,siz,£)Zp{z,E) 



fdE^^b,,s{z,£)Zp{z,E). (B.4) 



The first term vanishes since 6cas = for sufBciently low energies and Zp = beyond the maximal energy. 
The time integration of the l.h.s. between the present epoch (t = 0) and the first sources (t,nax) gives 



^max 

f d< [l.h.s.] = I dEEn^^{E)^. 



(B.5) 



hence we obtain Eq. ^. 



Appendix C. Comparison of Gamma- Ray Spectra 



Figure B.7 compares our calculations with spectra derived in Ref. [57l (Figs. 1 and 2) using n = 3 and 
7 = 2 (left plot) as well as n = and 7 — 2.6 (right plot), respectively, with (upper) and without (lower) 
BH contributions. 

The 7-ray spectra without BH pairs in the cascade are consistent with Ref. [57] apart from a slightly 
smaller energy density (recognizable as an overall shift downwards) probably due to the difference in the 
adopted IR/optical background. Since the energy density derived from our spectra agrees with the value 
obtained from Eq. ([6]) within 10%, we believe that the overall normalization of our spectra is correct. 

For a full calculation, i.e. including the BH pairs, our overall energy density agrees well with the calcu- 
lation in Ref. [57]. Again, the energy density derived from our spectra agrees with the value obtained from 
Eq. ([6| within 10%. The pair production dip at 10^ GeV is more pronounced in our spectra. Note that 
the energy loss of the cascade beyond lO'^ GeV is much more rapid than the Bethe-Heitler pair production 

rate, hence the modest increase in 7-rays beyond this energy meets our expectations. Moreover, instead of 

-7/4 



using a power-law approximation {An/<1E± ex E^ ) for the BH e spectrum (see Ref. [58j for a critical 
discussion) we use the exact differential cross-section of Ref. [51] . 
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Figure C.8: The effect of a large intergalactic magnetic field. The solid lines show the spectra of the best-fit for crossover 
energies -Emin = lO'^^'^ eV (left) and IQl^ eV (right) using B\q = 10~i^ G. The dashed lines show the corresponding results 
for a much larger field strength Big = 10~^ G. The 7-ray flux in the GeV-TeV region relevant for the Fermi-LAT spectrum is 
practically unaffected. 
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